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On beta distributed limits of iterated linear 
random functions 

Shaun McKinlay * 


Abstract 

We consider several special cases of iterations of random i.i.d. linear functions 
with beta distributed fixed points that generate nested interval schemes when 
iterated in a backward direction, and ergodic Markov chains in the forward 
direction. We prove that the fixed points are beta distributed by using a re¬ 
lated random equation with a gamma distributed solution that also generates 
a corresponding ergodic Markov chain with a gamma distributed stationary 
distribution. Our approach allows us to find limiting distributions of the 
random processes we consider, and provide solutions to the two random equa¬ 
tions, by solving only one of these equations. The paper extends many partial 
results available in the existing literature. 
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1 Introduction 

In this paper, we study two classes of discrete time stochastic processes that converge 
to beta distributed limits. The first class are interval-valued processes that are 
sequences of nested intervals lo = [0,1] D Ii D I 2 D • • •, where the length of the 
nth interval |I n | converges to zero almost surely (a.s.) as n —> 00 and the random 
rule of choosing interval I n+ i as part of is the same for each n > 0 (up to mirror 
reflections). The second class consists of time-homogeneous ergodic Markov chains 
taking values in [0,1]. 

The paper is centred around the following two key observations. Firstly, the two 
classes we consider are generated by the same family of linear random functions. A 
random function from this family generates an ergodic Markov chain (nested interval 
scheme) under forward (backward) iterations, which implies the following duality. 
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Every process in one class has a corresponding “dual” process in the other class, 
with the same limiting distribution (where the backward iterations converge a.s. 
and the forward iterations converge in distribution). This observation enables us to 
immediately obtain previously unknown limiting distributions of several processes of 
interest due to their duality with well studied processes appearing in the literature. 

Our second key observation is that one can show the limit is beta distributed by 
using a related random equation with a gamma distributed solution, allowing for 
simple analytic proofs. In addition, this random equation generates a corresponding 
ergodic Markov chain with a gamma distributed stationary distribution. 

These two key observations form the systematic approach used in the present 
paper, whereby we prove assertions concerning the limiting distributions of random 
processes we consider, and provide solutions to two random equations, by solving 
only one of them. 

Suppose Fi,F 2 , ... are i.i.d. random mappings of [0,1] into [0,1]. The forward 
iteration of the mappings is given by the compositions 

X n (-) := F n o F n _i o • • • o Fi(-), n > 1, X 0 (x) = x, (1) 

while the backward iteration of the mappings is given by 

Y n (-) := F 1 o F 2 o ■ ■ ■ o F n (-), n> 1, Y 0 (x) = x. (2) 

The theory of forward and backward iterations of i.i.d. random functions is well- 
established (see e.g. mum ED]). Here we will only summarise the results relevant to 
this study. 

It is clear that X n (-) has the same distribution as Y n (-) for each n. However, the 
properties of the forward and backward processes are very different. Under a broad 
contraction condition, for any fixed x E [0,1], the forward sequence {X n (a;)} n >o is an 
ergodic Markov chain, while the backward sequence {Y n (x )} n >o converges a.s. and is 
not in general a Markov chain (see e.g. [10]). In addition, the limiting distributions 
for the two sequences coincide due to the following contraction principle (see e.g. 
Proposition 1 in |4j). If Y := lim^,*, Y n (x) a.s. exists and does not depend on x (in 
which case we will say that iq is a contraction), then the Markov chain {X n (x)} n >o is 
ergodic with stationary distribution given by the law of Y. Following the terminology 
introduced in [2], we will refer to the forward and backward processes as being dual 
to each other. 

In this paper we are dealing with linear random mappings 

F n := .fA n .R ri i ra >1, (3) 

where fA n ,B n are random elements of the family 

& = ifa,b(x) = ax + 6(1 - x) : (a, b ) E [0, l] 2 }, (4) 

(H 1; Bi), (A 2 , B 2 ), • • • being an i.i.d. sequence of [0, l] 2 -valued random vectors with 
a given distribution //. 


2 


Iii that case, the forward process (CQ) is a special case of the random recurrence 

V n = D n V n - 1 + C n , n > 1, (5) 

where ( D n , C n ) are i.i.d. Revalued random vectors. If V n converges in distribution, 
the limiting distribution is the solution to the random equation 

V = DV + C, V independent of ( D,C ) = (Di,C\). (6) 

This is known as the perpetuity equation , since its solution may be interpreted as 
representing the present value of a commitment to pay at times n — 0,1,... random 
amounts C n+ 1 , subject to random discounting from period n to period n — 1 by 
D n (i.e. the present value of payment C n+ 1 is C\ for n = 0, and C n+ 1 n”=i Di 
otherwise). Due to this and numerous other applications (see for instance |5TI HU| ). 
the properties of recursion ([5]) have been studied extensively. A key reference is 
m, where sufficient conditions for the existence and uniqueness of the limiting 
distribution of V n as n —> oo are given. In particular, if the condition 

n 

^togl-Dfcl 4 -°° (7) 

k =1 

holds, then Lemma 1.5 in [40J implies that any solution V of (JHJ) is unique in dis¬ 
tribution and V n converges in distribution to V for all Vo- In the special case of the 
forward process (0Q) generated by the i.i.d. mappings (J3j) it is clear that condition 
O is satisfied when 

P(|Al — Bi\ = 1) < 1, (8) 

and this is precisely the condition required for F\ to be a contraction. 

Another notable paper is jl9j . where a multidimensional version of (J5J) is studied 
(i.e. C n and D n are d x d matrices, and V n is a d-dimensional vector). A number of 
studies have also considered special cases of (ED , where the stationary distribution 
is found by solving random equation (jUJ) (see mm and references therein), how¬ 
ever there are relatively few special cases where explicit formulae for the stationary 
distributions are known. The processes we consider below are generated by several 
special cases of perpetuity equation (J6D that admit closed form solutions that have 
not appeared in the literature. 

We will concentrate on the following special case of ()5]) that is generated by the 
forward iteration (JT]). Starting from a fixed x G [0,1] and writing X n = X n (x), the 
forward iteration is given by 

X n = (A n — B n )X n _i + B n , n > 1. (9) 

A deterministic version of the forward process (jUD (cf. ([I]), (13I) -(Hj)) was introduced 
by C. Li in the context of human genetics in f23j (where it was called the game of 
“give-and-take”). This process was subsequently studied (for non-degenerate /i) in 
[8], where conditions for convergence (in distribution) of (J9]) were given. In addition, 
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[8] included several examples of distributions p on [0, l] 2 generating Markov chains 
by dnj) with known stationary distributions. More recently, a multivariate version of 
this give-and-take model was considered in [25] . 

If Fi is a contraction, the unique stationary distribution P of the forward process 
(jSJ) satisfies the random equation 

X = AX + B( 1 - X), X ~ P independent of (A, B) = (A 1; B{), (10) 

where A" ~ P denotes that X has distribution P. Therefore, the stationary distri¬ 
bution of the forward process and the limiting distribution of the backward process 
coincide with the law of solution to (HU . 

We can interpret the forward process ([T]) in terms of the movement of a particle in 
[0,1], where X n is the location of the particle at time n, and the particle moves from 
X n _i to X n by (jH]). In particular, we will consider models where for each n either 
A n = 1, corresponding to a move to the right, or B n = 0, corresponding to a move to 
the left. Several special cases of models of this type (where, at each step, a particle 
at X n G [0,1] moves at time n + 1 either left to a random point in [0, A n ], or right 
to a random point in \X n , 1]) have been shown to have beta distributed stationary 
distributions (see e.g. [8l uni EH [38]). See also [321 [26] . where an extension of this 
model to the case when the direction of the next move is a function of the particles 
current location was considered. 

We will interpret the backwards process (J2[) as generating a sequence of nested 
intervals lo = [0,1] D Ii D I 2 D ■ ■ ■ given by the ranges of the corresponding 
backward mappings lo, Ti, I 2 , • • • (he. Y n maps [0,1] onto I n ), and consider only the 
case when |I n | — y 0 as n — )■ 00 (i.e. when F\ is a contraction). A nested interval 
scheme of this type was studied in nang, where the limiting location Y was shown 
to be beta distributed. 

In some of the cases we consider, the nested interval scheme generated by the 
backward iteration ([2]) is an interval splitting scheme of the following type. Choose a 
random (not necessarily uniformly distributed) splitting point Si in [0,1] and select 
according to a given (possibly random) rule one of the subintervals [0, Si], [Si,l]. 
Continuing this procedure in the same way independently on the chosen subinterval, 
we obtain a random sequence {S n } n >i. When the law of the splitting point Si is not 
concentrated on {0,1}, the length of the chosen interval tends to zero a.s. as n —> 0, 
and therefore S n converges a.s. to a random variable Y. In several cases considered 
previously, this limiting random variable Y turns out to be beta distributed (see 
e.g. [SI 0 0 [El CEE]). We call this scheme the CGZ model after Chen, Goodman 
and Zame who first studied it in [6] under the conditions that S n ~ I/[0,1] and the 
largest subinterval is selected with probability p at each step. 

Our systematic approach allows us to describe the nested interval scheme and 
ergodic Markov chain that are dual to each other in each case we consider below. 
In particular, we will describe the dual of the above process when the larger of 
the subintervals [0, Si], [Si, 1] is chosen with probability p, and the smaller with 
probability 1 — p, that eluded the authors of [2] . 
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Most of the nested interval schemes considered in the literature can be classified 
using a convenient notation of the form 



where L , LL", M 7 , M" are letters and N is a natural number which have the fol¬ 
lowing meaning. 

Imagine each step of a nested interval scheme I 0 = [0,1] D Ii D I 2 D ■ ■ ■ is 
generated by the following two stage process. In the first stage of step n we choose 
a fixed (one and the same for all n) number of subintervals of I n _i, n > 1. In the 
second stage, we select one of the chosen subintervals to be I n . 

We set L := S if the intervals from the first stage can be constructed by splitting 
the interval I n _i, n > 1 , into two or more subintervals (i.e. they form a partition of 
I n _i), otherwise we set L := G for general. The number N is equal to the number 
of intervals generated in the first stage. We set L'D if the rule for the first stage 
is deterministic , or L := R if the rule is random. The subscript M' of the second 
component is equal to / if the rule is the same (i.i.d. after scaling and translation) 
for each stage of the process, otherwise M' := G for general. Similarly, in the final 
component we set L" D (resp. R) if the rule for the second stage is deterministic 
(resp. random), and M" := I if the rule is the same (i.i.d.) for each stage of the 
process, or M" := G otherwise. 

We will mostly be interested in nested interval schemes with M' = M" = / since 
in that case the scheme is completely determined by the rule for choosing the first 
interval Ii from [0,1] (e.g. the CLZ model is an (S 2 , Ri\Ri) scheme in the general 
case). 

To illustrate our classification, consider the following scheme that was introduced 
by Kennedy in the elegant note id- Suppose at the nth stage of the process 
I n = [A n , B n ]. Take k i.i.d. points u[ n \ ..., uniformly distributed on [A n , B n ], 
and let C n := , D n := , where < ■ ■■ < are the order statistics 

of u[ n \ ..., U^\ Then the (n + l)st interval I n+ i is taken to be [C n , B n \, [A n , D n \ 
or [C n , D n \ independently at each stage with probabilities p,q,r, respectively, with 
p + q + r = 1. Then, as shown in na one has A n , B n Z ~ f3(k(p + r), k(q + r)). 

Consider the case p,q,r > 0. For k > 3 the above scheme is an (G 3 , Ri\Ri) 
nested interval scheme, whereas for k = 3 it is an (S 3 , Rj\Rj) scheme. Note that 
interval schemes can admit different classifications within the introduced notation. 
For example, when p — 1, Kennedy’s scheme can be viewed as either (S 3 , Rj\Dj) 
or (G\, Rj\Dj) , where in the latter we only considered the middle interval [C n ,D n ], 
rather than choosing the middle interval deterministically from the three intervals 
in the former. Despite that, the classification may still be useful for the schemes we 
will consider. 

Further, we observe that the forward and backward process are also related to 
products of i.i.d. 2 x 2 random stochastic matrices since the processes can be written 
equivalently by 



Ax 1-A 

B 1 1 -B 1 



A n 1 — A n 

Bn 1 -B n 


) . ( 11 ) 
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and 


(y.(*),i-y.(*)) = (*,i-*)(^ LbI)'"(b; !-£)■ (12) 

respectively. If F\ is a contraction, then it is clear that for the left product we have 
that 

ft !:£)■•■(£ 1:10 < 13 > 

as n —)■ oo, where Y := linp^oo Y n (x) a.s., and for the right product we have that 
(A 1 1 - AA (A n 1 -A n \ d (X l~X\d z (Y 1 - Y\ 

V^i 1 1-bJ {x 1-xJ [y 1-y) 

as n —y oo, where X satisfies the random equation (ITU]) that can be equivalently 
written as 

(V 1 - V) = (V 1 - V) (g Jig). (14) 

It is well-known that a necessary and sufficient condition for the convergence (1T51) 
is that n is not concentrated on {(1, 0), (0,1)} (see [13]), which is equivalent to our 
contraction condition ([51). 

Rather surprisingly, in the literature devoted to the study of Markov chains 
generated by © and nested interval schemes generated by (EJ), the authors mostly 
do not note and/or use the fact that the forward and backward iterations can be 
represented by the right and left products of i.i.d. 2x2 random stochastic matrices, 
respectively. One exception is [27], where rates of convergence for the left products 
(1151) are given. Observing that the nth stage of the interval splitting schemes that 
were considered in isunEun] correspond to the left hand side of (IT5|) . the author 
of [27] gives geometric convergence rates in terms the minimal L p -metric and also 
(where possible) for the Kolmogorov (uniform) metric. We note that the same 
approach may be used to derive similar convergence results for the special cases we 
present below. 

Interestingly, simply recognising that the stationary distribution of the forward 
process (P is the law of the almost sure limit in (fT5j) would have progressed the 
search for special cases of the latter that admit closed form expressions for those 
distributions. To illustrate this observation denote by (3(a,b), a, b > 0, the beta 
distribution with density 

B(a, b)~ 1 u a ~ 1 ( 1 — u) b ~ l , 0 < u < 1, 

where B(a , b ) is the beta function. In addition, if A ~ Fa is independent of B ~ F B , 
we write Fa <E> Fb for the law of ( A, B ). Then, one can observe that the main 
assertion (Theorem 2) in [39| devoted to random matrix products assumed that 
H = (3(a, a) <S> (3(a, a), in which case Y ~ (3(2a, 2a) for the limit on the right hand 
side of (1151) . However, more than two decades earlier in paper [8] on random linear 
iterations, it was shown that for /i = (3(a, b) ® (3(b,a ), one has that X n (x) 4 I ~ 
(3(a T b, a T b). 
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The remainder of the paper is structured as follows. In Section [2j we present our 
main results, beginning in Section I2TT1 with a description of the general nested interval 
scheme generated by the backward iteration ((2]) . We then devote Section 12.21 to 
proving an equivalence result between beta distributed solutions to (|IU]) and gamma 
distributed solutions to a related random equation. In Section 12.31 we introduce 
several equivalent claims concerning limiting distributions for processes and solutions 
to random equations we consider. 

In Sections 12.41 and 12.51 we use the approach discussed above to find the explicit 
form for the limiting distributions of nested interval schemes and Markov chains 
on [0,1] in two broad special cases. In the first one, at each step either A n = 1 
or B n = 0. In the second more general case, the distribution of (A n , B n ) can have 
support on an arbitrary subset of the unit square [0, l] 2 , with most cases considered 
being closely related to the beta distribution. 


2 Main Results 

We start by describing the general nested interval scheme generated by the backward 
iteration ([2]). 


2.1 A general nested interval scheme 

Recall that F\. F 2l ... are the random functions from & generated by a sequence 
i.i.d. vectors (Ai,Bi), (A 2 ,B 2 ),... with distribution p (see (J3J), ([I])). Plainly, 

Fi(x) = (Ai — Bi)x + Bi 

maps [0,1] onto Ii = [min(d 1 , Bi), max(d 1 , iR)]. So the nth backward iteration 
from (J2D maps [0,1] onto the nth interval of the following nested interval scheme. 
Starting with the unit interval [0,1], select the interval lb with end points A\ and 
B\. For the second step, noting that the mapping F 2 is linear, we see that the end 
points of the second nested interval I 2 C Ii are given by 

A 2 := 1^(1) = A 2 (Ai — Bi) + B\ and B 2 : = W(0) = B 2 (Ai — B\) + B\. 


It follows that 


/ {A 2 — B 1 , B 2 — F?i) 

V A 1 -B 1 

~ ((A 2 — A, B 2 — A\) 

V B l -A l 


A 1 -B 1 >0 
A - B 1 < 0 


£((A,£i)), 


£((1 -A 1 ,l-B 1 )), 


(15) 


where jC(X) denotes the distribution of X. This is, if A\ > B\, the law of the 
ends points of the interval I 2 conditioned on lb is the same (after translating and 
scaling to the unit interval) as that of the end points of lb. Otherwise, if Ai < Bi, 
the law of (^ 2 ,^ 2 ) conditioned on Ii is the same (after translating and scaling) as 
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that of (1 — Ai, 1 — B i). Similarly, the end points of the nth interval are given by 
A n \= Y n ( 1) and B n := Y„(0), with 



£((A,5i)), 


(16) 




We will mostly be interested in special cases where one of the following conditions 
is met: 


[Ml] A > B a.s., 

[M2] // = //':= £((1 - A, 1 - B)), 

[M3] M ' = £((B,A)). 

In that case the nested interval schemes generated by ([2]) satisfy the following i.i.d. 
property. 

[I] For all n > 1, the conditional distribution of 


In - min I n _i 


given I n _i doesn’t depend on I n _i and is one and the same for all n > 1. 

2.2 A relationship between beta and gamma distributed so¬ 
lutions to linear random equations 

The following proposition is fundamental to our approach and provides a simple 
way to show that solutions to (TTUT) are beta distributed. Denote by T a the gamma 
distribution with scale parameter 1 and shape parameter a G M + := (0,oo), with 
density n a_1 e _ “r(a) _1 , u > 0. 

Proposition 1. Random equation (fTUj) has solution X ~ /3(a,b) for some a,b > 0 



holds for {V \, V 2 ) ~ T a <8) r 6 . 


Proof. The proof closely follows that of Theorem 3 in [25]. Suppose that relation 
(1T7D holds. Then 


AV t + BV 2 
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Pi + P 2 


P1 + P2 


(18) 







where V ~ T a+ b is independent of ( A , B , hi, V 2 ), the equality in distribution follow¬ 
ing from Lemma 1 in [25]. Similarly, we have that 


V, = 


V! 


Vj+Vj 


(W + V 2 ) = 


V! 


V, +v 2 


V, 


(19) 


where V ~ r a+ h is independent of {V\ , V 2 ). Taking logarithms of the right-hand 
sides of (TTST) and f[T9|L we obtain 


In 


(AV X +BV 2 \ 

V V 1 + V 2 J 




( 20 ) 


Denoting by -0, ip and x the characteristic functions of the first, second (and fourth), 
and third terms in (I 2 U|) . respectively, we have 


^(u)<p(u) = x(u)<fi(u)- 

It is not hard to show that (p(u) 7 ^ 0. Therefore ij) = %, and so 

hi d AV 1 + BV 2 
hi + V 2 - hi + V 2 ■ 


We conclude that random equation (flUjl has solution Y = hi/(hi + V 2 ) ~ /3(a,b). 

Conversely, suppose that random equation (fTUj) has solution Y ~ /3 (a, b). Follow¬ 
ing the same steps as above in reverse order, we conclude that (fT7j) holds. Proposition 
Q] is proved. □ 

Proposition Q] relates the original sequences of iterations of functions on [0,1] 
to the following forward and backward iterations defined on R + . For the family 
{g a ,c( x ') = ax '+ c: i a , c ) £ [ 0 , 1] xM + } and an i.i.d. sequence (A 1; Bi, hi), (A 2 , B 2 , hi), 
with distribution p <S> r&, set G n := gA n ,B n v„ ■ Denote the forward iterations of 
{G n } n > 1 by 


X' n (-):=G n oG n _ 1 o...oG 1 (-), n> 1, X' 0 (x)=x', (21) 

and its backward iterations by 

r/(-) :=G 1 oG 2 o...oG n (-), n> 1, Y'(x) = x’. (22) 

Starting from a hxed x' > 0 and writing X' n = X' n (x'), the forward process (l2l| is 
given by 

X'n ■= A n X' n _ 1 + B n V n , n > 1. (23) 

This is also a special case of recursion (|S]), and so X' n converges in distribution to 
a unique (in distribution) X' when condition (JTJ) is satisfied. In our case, since 
Ai G [ 0 , 1 ], this conditions is satisfied when 

P(kLi = 1) < 1. (24) 
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2.3 Some equivalent claims 

For (a, b) 6 Rl, we introduce the following claims to be used in the statements of 
assertions in the remainder of this paper: 

[All Suppose V 2 r 5 . Then V\ ^ T a is the unique solution to (1 1 7f) . 

[All One has X ~ j3(a, b ) is the unique solution to (fTUjb 

[A2] For Y n {x) given by (j2J), as n —> oo, one has Y n (x) -^4- Y ~ /3(a,b ) uniformly 
in x G [0,1]. 

[A3] For X n (x) being the location of the particle after the nth step of the Markov 
chain ([I]) starting from x G [0,1], one has X n (x) -4 X ~ /3(a, b) as n — * oo. 

[A4] For X' n (x') being the location of the particle after the nth step of the Markov 
chain (12T|) starting from x' G R+, one has X' n (x') -4 X' ~ T a as n —> oo. 

The next proposition shows that to establish [Al']-[A4], it suffices to just prove 
[Al 7 ] (or [Al 77 ]). 

Proposition 2. Suppose that the distribution p of (Ai, Bi) is such that the conver¬ 
gence conditions (jHJ), <l24|i hold, and at least one of [Al 7 ] or [Al 77 ] holds true. Then 
all assertions [A1']-[A4] hold true. 

Proof. Since the convergence conditions (jHJ) , (l24j) are satisfied, any solutions to ran¬ 
dom equations ffTOj) and (ITT)) is unique. Therefore by Proposition [T| [Al 7 ] holds if 
and only if [Al 77 ] holds. Further, the law of V\ from [Al 7 ] is the stationary distribu¬ 
tion of the Markov chain from [A4], and the law of X from [Al 77 ] is the stationary 
distribution of the Markov chain from [A3] which coincides with the law of the lim¬ 
iting location in the nested interval scheme from [A2] by the contraction principle 
in [4] . □ 

We will now consider several classes of distributions p on [0, l ] 2 satisfying con¬ 
vergence conditions (jS]), (ITT)) , with beta distributed solutions to random equation 
ca. Uniqueness of the solutions to ca and (fTTj) is implied by the convergence 
conditions. 

2.4 Model 1 

In this model we assume that P({A = 1, B < 1} U {A > 0, B = 0}) = 1. Then for 
some fixed p G (0,1), one can write 

(A, B) 4 4(1 - L, 0) + (1 - 4)( 1, R ), (25) 

where I p ~ Be(p ) is independent of the random vector (L, R ) G [0, l ) 2 and Befp) 
denotes the Bernoulli distribution with success probability p. 
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This model has two important features. Firstly, since A > B a.s., the nested 
interval scheme generated by (J2J) satisfies the independence property [I]. Secondly, 
the representation (I25j) allows for a simple description of the nested interval scheme 
generated by (]2j), the Markov chain on [0,1] given by (IHj) , and the Markov chain 
on M + given by Q23j) . These processes are described respectively in (1A), (IB), and 
(1C) below. In particular, the nested interval scheme ((G2, Ri\Ri) in the general 
case) described by (1A) and the Markov chain described by (IB) are dual in the 
sense of the contraction principle (Proposition 1) in [4]. 

(IA) Starting with the unit interval, we choose the two subintervals [0,1 — L\ and 
[. R , 1], i.e. by moving the right (resp. left) end point of the interval to the left 
(resp. right) a distance given by the proportion L (resp. R) of the interval. 
Next, select the interval [0,1 — L] with probability p, or [ R , 1] with probability 
1 — p, continuing independently in the same way on the chosen subinterval. If 
1 — L = R, then this is the CGZ model ((£ 2 , Ri\Ri) in the general case). 

(IB) A particle starts from a point x E [0,1] and moves left the distance Lx with 
probability p, or right the distance R( 1 — x) with probability 1 —p. The process 
continues in this way, each new step being independent of the past. 

(IC) A particle starts from a point x' E M+ and moves left the distance Lx' with 
probability p, or right the distance RV\ with probability 1 — p, where V± ~ T a 
for some a > 0. The process continues in this way, each new step being 
independent of the past. 

Clearly assertions [A2]-[A4] give the limiting distributions of the processes de¬ 
scribed in (1A)-(1C), respectively. 

The first special case of Model 1 considered below demonstrates two points. 
First, it is not immediate that the backward and forward process are related when 
studying either independently. Indeed, a search through the literature found only 
one previous study (see m) that examines both the nested interval (backward) 
process and the ergodic Markov chain (forward) process generated by distributions 
/i on [0, l] 2 considered in the present paper. However, [TOj only discussed the general 
properties of the backward process (i.e. a.s. convergence of Y n (x) for all x E [0,1]). 
Second, it is often possible to gain additional insights into the limiting distribution 
by considering both the forward and backward process together, rather than only 
one in isolation. 

Consider the following (S 2 , Dj\R/) version of the CGZ model. Fix a p E (0,1) 
and, starting with the unit interval, choose the interval [0, p] with probability p, 
or [p, 1] with probability 1 — p. The next nested interval is obtained in the same 
way from the first chosen interval, and so on. It is clear that this scheme can be 
generated by a single uniform random variable U ~ U[ 0,1] as follows. On the first 
step, of the two intervals [0,p] and (p, 1], we choose the one containing the value of 
U. At the next step, we partition the chosen interval into two with the same ratio 
of lengths, and again choose the subinterval containing the value of U, and so on. It 
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is obvious from this representation that the limiting location of the nested intervals 
has distribution {7[0,1]. 

Thus [Al"] holds for this scheme. It follows from Proposition [2] that [Al'] will 
also hold. That fact however admits an elementary direct proof as well. 

Theorem 1. Under assumption f|25p . let L = (1 — p), R = p. Then [Al'] holds for 
a = b = 1. 

Proof. Since the V)’s are independent exponential when a = h = 1, one has 
Ee it(AV 1 + BV 2 ) = pEe itpV 1 + (1 _ p)Ee im+pV2) 

= p( 1 — itp)~ l + (1 — p)(l — i£) -1 (l — itp)- 1 = (1 — it) -1 , 

which is the characteristic function of V± as required. □ 

Note that assertion [A3] under the conditions of Theorem [T] extends Theorem 2 
in [3Sj where it was proved in the special case p = 1/2. Likewise, assertion [A2] in 
this case when interpreted for the products of 2 x 2 stochastic matrices (cf. ffl3]) ) 
extends the result of example on p.160 of [33] that also considered the special case 
when p — 1/2. 

The next case extends the situation considered in [37] where assertion [A3] from 
Theorem [7] below was first proved, extending that assertion in the case z — 1 from 

m- 

Theorem 2. Under assumption (|25]h let L = R ~ (3(1, z) for a fixed z > 0. Then 
[A1']-[A4] holds for ( a,b ) = ((1 —p)z,pz). 

Proof. According to Proposition (]2]) it suffices to show that [Al"] holds. Clearly 
random equation (ITO can be written as 

X = (1 - C)X + CD, (26) 

where X is independent of 

(C, D) = (1 - A + B, B/( 1 - A + B)) =| 1^°/’ / I g’ (27) 

and (C, D ) has distribution 

(C,D) ~P(l,z) ®Be(l -p). (28) 

Therefore, it suffices to show that X ~ /3((1 — p)z,pz) is the unique solution to 
(l26j) with ( C, D) distributed as per (l28jh Uniqueness follows from Lemma 1.5 in [40] 
since condition ((Jj) is satisfied for an i.i.d. sequence with same law as (1 — C). The 
fact that X ~ (3(( 1 — p)z,pz) is the solution follows from Section 3 in [32] (see also 
Proposition 1 in [37] where it was shown that (125]) holds for X £((1 ~P)z,Pz) 
using the method of moments). □ 
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Remark 1. Observe that for z — 1, assertion [A2] under the conditions of Theorem 
[2] gives the limit distribution of the following (S 2 , Rj\Rj) version of the CGZ model 
that was considered in [IB]. Choose a uniformly distributed point U\ in [0,1], and 
select (independently of Ui) the subinterval [ 0 , Ui\ with probability p, or [f/i, 1 ] 
with probability 1 — p. Continuing independently in the same way on the chosen 
subinterval, the limiting location of this interval splitting scheme has distribution 
fj{l —p,p). This result is an extension of the main result in [7] where it is proved 
for p — 1/2. The authors of [7] describe the fact that the limit distribution of that 
scheme is the arcsine law as “very startling and amazing”. The next remark shows 
that the result obtained in [7] follows from arcsine laws of the Brownian motion 
process and the uniform laws of the Brownian bridge process. 

Remark 2. For a standard Brownian motion W = { W t : t > 0} started at zero, 
Levy’s “First Arcsine Law” (see e.g. p.136 in [28|) asserts that the time 

T\ := inf {t : W t = max W s \ 

0<s<l 

of the first maximum on [ 0 , 1 ], and the time 

r 2 := sup{t G [0,1] : W t = 0} 

of the last zero of W on [0,1], are both arcsine distributed. Levy’s “Second Arcsine 
Law” (see e.g. p.137 in [28j) says that the same arcsine law holds for the occupation 
time of the positive half axis of W on [0,1]. That is, for 

u(t) := / 1 { m 4 >o }ds, 

Jo 

where 1 a denotes the indicator function for event A, one has z/(l) | 8 ( 1 / 2 , 1 / 2 ). 

We will now briefly summarise the approach used by P. Levy in [22] to prove 
the Second Arcsine Law (for a more detailed outline of this approach, see e.g. [29]). 
First, Levy showed that, given r 2 = s, the process up to time s (i.e. { W u : 0 < 
u < s}) is a Brownian bridge of length s. Observing that the time spent on the 
positive half axis by a Brownian bridge of length s is uniformly distributed on [0, s] 
(in fact, this property holds for a much broader class of bridge processes having 
cyclically exchangeable increments, cf. 0), he derived the arcsine law for u(l) using 
the uniformity of v(s)/s and the fact that after s the Brownian path is equally likely 
to be positive of negative. 

That is, Levy showed that t/r 2 + /i/ 2 (l — t 2 ) is arcsine distributed, where 
(17,/i/ 2 ,r 2 ) ~ f/[0,1] ® Be( 1/2) (8) /3( 1/2,1/2), and this is precisely assertion [Al"] 
under the conditions of Theorem [2] in the special case (z,p) = (1,1/2). 

Note that the same approach can be used to show that {y{l)\Wi < 0) ~ 
/3( 1 / 2 , 3/2) (and similarly that > 0) ~ /3(3/2, 1 / 2 )). Arguing as above, it 

suffices to show that 

Gt 2 ~/3( 1/2, 3/2) 
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for (U,r 2 ) ~ t/[0,1] ® (3( 1/2,1/2), which follows from the well-known fact that 
AB = C for ( A , B) ~ (3(a + b, c) ® /3(a, b ), and C ~ (3(a, b + c). 

Remark 3. It was suggested in [25] that “maybe the only way” to prove the special 
case of assertion [A3] under the conditions of Theorem [2] when z — 1 (cf. Theorem 
3 in [38]) is via the moment convergence theorem. Here we will give another proof 
Theorem [2] using a result from the well-known Polya urn model. 

In the that scheme, a single urn contains w white balls and b black balls, and 
at each stage a ball is drawn at random, and then returned to the urn with one 
additional ball of the same colour. If there are b n black balls and w n white balls in 
the urn after n draws, then for the proportion P n := b n /(b n + w n ) of black balls in 
the urn, one has P n+ k P% r-u (3(b n ,w n ) (see e.g. p.243 in [13]). 

Suppose the urn contains w white balls and b black balls. Since after the first 
draw a white (black) ball is chosen with probability w/(w + b) (b/(w + b)) it follows 
from the above observation about almost sure convergence that Pn +1 ~ W + 1) 
(p :+i (3(b + l,iu)). We therefore have that 

x = I b/{b+w) x b + (1 - I b/{b+w) )X w , (29) 

where A" ~ (3(b,w) and I b /( b + w ) ~ Be(b/(b + w)) is independent of X b ~ (3(b + l,iu) 
and X w ~ /3(b, w + 1). Now, suppose that (Vi, 14, V3) ~ Ti <E) T^ <g) T^, then (129|1 can 
be written equivalently as 

^2 d ^1 + ^2 ( , V 2 

v 2 + 1/3 b/{b+w) v 1 + v 2 + v 3 + { b/{b+w) ) v 1 + v 2 + v 3 

v 2 v 2 + h 3 ki 

V 2 + Hi + V 2 + 14 b/[b+w) Vi + V 2 + V 3 

d v 2 ^2 + ^3 v; 

~V 2 + V 3 v 1 + V 2 + V 3 b/{b+w) Vl + Vi + V^ 

where the last equality follows from Lemma 1 in [25], and (Vi, V 2 , V 3 ) is an indepen¬ 
dent copy of (Vi,V 2 ,V 3 ). The right hand side of the last line above has the same 
form as the right hand side of (125]) for (C, D) ~ /3(1, b + w) <g) Be(b/(b + w)), with 
solution X P(b,w). 

One can also show that Theorem 2 holds by proving assertion [Al'] analytically 
using the method of characteristic functions. I11 doing so now, we will establish the 
limiting distribution in two cases that have not been considered previously in the 
literature from the well-known properties of the Gauss hypergeometric function 

2 Fi(a,b] c\rj) = B(b,c — 6)” 1 f u b ~ l (l — u) c ~ b ~ l (l — ur])~ a du, 

Jo 

for 77 6 C, Re(c) > Re(fe) > 0, and |arg(l — 77)| < n (see e.g. 15.3.1 in [T]). Now 
for (14,14) ~ T(i_ p )^ <8) T p „, L = R~ (3(1, v) and (A, B) given by (125]) . and letting 
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z — (1 — p)v , y = pv , we have 
Ee it(AV 1 +BV 2 ) 

= pEU t{l ~ L)Vl + (1 - p) Ee it{Vl+RV2) 


y 


y + z 
z 


+ 


y + z 

y i 


y + z 

z 


E ( e i t(i-i)v 1 |( 1 _ = U )P((! _ L) G du) 

3 

Ee i4Vl / E(e™|R = u)P(i2 G du) 

Jo 

C (1 - itu)~ z u y+z ~ 1 B(y + z, l)' 1 du 


H--—(1 — it) 2 / (1 — itu) y (l — u) y+z 1 B(l,y + z) 1 du 


y + z 

y 


y + z 

y 


2 F l (z,y + z-y + z + l-,it)-{ --—(1 - it) 2 2 Fi(y, 1; y + z + 1; it) (30) 


y + z 


2 Fi(z, y + z\y + z + l\it) + 


y+z y+z 

= 2 Fi{z, y + z + 1; y + z + 1; it) 

= (! - 


2 F 1 (z + l,y + z;y + z + l]it) (31) 

(32) 

(33) 

where (j3lil . f[32]l . and (1331) follow from 15.3.3, 15.2.14, and 15.1.8 in [I], respectively. 
The right hand side of (l33l) is the characteristic function of Vi, and so [Al 7 ] holds 
as claimed. 

We will now use the expression (l30l) and the well-known property that 


2 Fi(a,b;c;rj) = 2 Fi{b, a; c; rj) 


(34) 


to establish the limiting distribution of our next case. Indeed, using this property 
it is clear that expression (l30]l is equal to 

—7— 2 Fi(z,y + z-,y + z + l;it) H- J— (1 - it)~ z 2 F 1 (1, y, y + z + 1 ;it), (35) 

y+z y+z 

and this implies that assertion [Al'] holds under assumption (I25|l with L /3(l ,y + 
z), R ~ j3(y, z + 1), p — y/(y + z), and (V \, V 2 ) ~ T y <8) T^ We have proved the 
following theorem. 

Theorem 3. Under assumption (125]) . let L ~ /3(l,y + z), R ~ fd{y,z + 1), and 
P — y/(y + z )> y-, z > o. Then [A1']-[A4] hold for (a, b ) = (z, 1). 

We can derive yet another case from (135|) as follows. Recalling that (133|) is equal 
to (135]) . we have 

(1 -ity z = — ^{z^y + z\y + z + l;it) H-(1 - it)~ z 2 T\(1, y, y + z + 1; it) 

y+ z y+z 

= -^(l-it) 1 ~ z 2 F 1 (l,y + l;y + z + l]it) 

y + z 

H-— (1 - it)~ z 2 Fi( 1, y, y + z + 1; it), 

y + z 
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the second equality following from relations (l34|i and 15.3.3 in [1), and therefore 

(1 - ity L = —— 2 ^( 1 , y + 1 ; y + z + 1 ; it) + —(1 - it ) _1 2 Fi ( 1 , y , y + z + 1 ; it), 
y+z y+z 

The last equality implies that assertion [A1 '] holds under assumption (l25]i with 
L ~ /3(z,y + 1), R ~ P(y,z + 1), p = y/{y + z), and (Vi,V 2 ) ~ T 1 ® Ti. The 
following theorem is proved. 

Theorem 4. Under assumption (l25j) . /ef L P(z,y + !), ^ /?(?/, z + 1), and 

P = y/(y + z )> y-i z > o. Then [AT]-[A4] hold for a = b = 1. 


Remark 4. Note that assertion [A3] under the conditions of Theorem 0 is an 
extension of Proposition 2 in [37], where it was proved in the special case y = z. 
In that case, we can derive an interesting distributional relation involving mixtures 
of gamma distributions from which we will obtain an apparently new explicit form 
solution to the perpetuity equation as a special case. 

Indeed, using the equivalent random equation (1251) we have from assertion [Al"] 
under conditions of Theorem 0] that 

U = U{l-C) + h /2 C , 


for (17, C, I 1 / 2 ) ~ U[ 0,1 ]®(3(y, y+l)®Be(l/2). Therefore, for (Vi, V 2 ) ~ Y y / 2 ®T y / 2+1 


U{ V, + V 2 ) = 

d 


UV 2 + h /2 v x 
Vh + K, 
UV 2 + I X/2 V X 

V\ T V2 


(V! + V 2 ) 
(V! + V 2 ) 


UV 2 + Iy 2^1, 


(36) 


where (Vi, V2) is an independent copy of (Vi, V2), and the second equality in distri¬ 
bution follows from Lemma 1 in [2.5]• 

Let {7(t)}t>o, 7(0) = 0 be a gamma process with independent increments such 
that 7 (t) — 7(5) ~ T 4 _ s for 0 < s < t. Suppose also that {71 (-)}? { 72 ( - )}> • • • are 
independent copies of the process {7(t)}t>o, U ~ U[0, 1], and Ji, J 2 ,... are i.i.d. 
with distribution Re(l/2). Assuming independence of all the processes and random 
variables introduced above, we can rewrite random equation (l3fi]i as 


U^(y + 1) = U^(y/2 + 1) + J^y/2), 


and, by repeated substitution, we obtain 

OO 

U 7 (y + 1) = f/ 7 (l) + JkUk(y/2 k ), 

k=\ 

where the series converges a.s. by Kolmogorov’s two series theorem (see e.g. p.386 
in [35]). By conditioning on {Jk} we see that 

OO / OO j \ 

J klk{y/ 2 k ) = 71 ( y -£ ) = 'YiiyUf), 

k= 1 V k= 1 Z J 
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where U\ ~ £7[0,1] is independent of ({7(-)}, (7i(-)}, U). So we obtained the curious 
relation 

U'y(y + l)^U'y(l)+'y 1 (yU l ). (37) 

In particular, for y — 1, the random variable 

f/ 7 (l)+ 7l (f/ 1 ) 

is exponentially distributed, and therefore V ~ U[ 0,1] is the unique solution to the 
perpetuity equation ([6]) with (D,C) = (U, 7 i(£/i)). 

One can also show that (137[) holds by computing the Laplace transform (or 
characteristic function) of each side. A search of the literature yielded only one 
paper [IB] that considered the mixture , yi(yUi). In that paper the explicit form for 
the density of y(I7)/7(l) was given. 

The well-studied original (S 2 , Ri\Ri) version of CGZ interval splitting model 
from [6] is as follows. Choose a uniformly distributed point U in [0,1], and select 
the larger of the subintervals [0, U], [U, 1] with probability p, and the smaller with 
probability 1 — p, and continue this procedure independently on the chosen interval 
and so on, denoting the limiting location by Z p . In that paper the authors show 
that Z\ ~ /3(2, 2), and derive several properties for the distribution of Z 0 . The 
related paper [7] is devoted to proving that Zy 2 ~ /3( 1/2,1/2) (note that this 
result is a special case of the assertion [A2] under the conditions of Theorem [2] for 
(z,p) = (1,1/2)). These results were obtained independently in the later paper 
P], that also showed that Z p is beta distributed if and only if p e {1/2,1}. More 
recently, the problem of finding the dual of this scheme motivated the authors of [2] 
to consider three new schemes, neither of these being the dual. 

It is easy to verify that the scheme above is generated by the backward iteration 
([2]) with F n (x) of the following form: 

1{ u>i/ 2 }[I P Ux + (1 - I p )(x + £7(1 - x))] + l{c/<i/2 }[I P (x + U{ 1 - x)) + (1 - I P )Ux], 

where (I/, I p ) ~ U[0, 1] ® Be(p). We will now describe the dual process of the above 
CGZ model, that eluded the authors of [2]. 

A particle starts from a point x G [0,1] and moves left or right with equal 
probabilities. If the particle moves left, it moves either to a uniformly distributed 
point in [x/2,x] with probability p, or to a uniformly distributed point in [0,x/2] 
with probability 1 — p. Similarly, if the particle moves right, it moves either to a 
uniformly distributed point in [x, x+ (1 — x)/2] with probability p, or to a uniformly 
distributed point in [x + (1 — x)/2, 1] with probability 1 — p. The process continues 
in this way, each new step being independent of the past. 

It is clear that 0, 1/2 and 1 are the only values of p resulting in the particle 
moving to a point uniformly distributed on a subset of [0,1]. The three cases are 
described below: 

• p — 0: The particle moves to a uniformly distributed point on the furthest 
half of the interval on the left or right of x (i.e. the intervals [0,x/2] or [x + 
(1 — x)/2, 1], respectively). 
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• p = 1/2: The particle moves to a uniformly distributed point on the entire 
subinterval of [0,1] on its left or right. 

• p — 1: The particle moves a uniformly distributed point on the closest half of 
the interval on the left or right of x (i.e. to the intervals [x/2, x\ or [x, x + (1 — 
x)/2], respectively). 


Recall that the limit is only beta distributed when p £ {1/2,1} (cf. [Oj), and 
observe that for p = 1/2 the process corresponds to the case where p' = 1/2 for the 
following scheme. A particle starts from a fixed x £ [0,1] and moves left (resp. right) 
with probability p' (resp. 1 — p') a distance given by a (3( 1, z) distributed proportion 
of the length of [0, x\ (resp. [x, 1]). The particle continues in same way independently 
of the past, with limiting distribution /?((1 — p')z,p'z) (cf. assertion [A3] under the 
conditions of Theorem [21) - It is natural to ask if the limit also beta distributed when 
the particle moves a distance given by a {3(1, z) distributed proportion of the length 
of the closest half subintervals to the left or right. It happens that if p' = 1/2, then 
the limit is indeed beta distributed, as we will show in the following theorem. 

Theorem 5. Under assumption (1251) . let 2 L = 2 R ~ {3(1, z), z > 0, and p = 1/2. 
Then [Al']-[A4] hold for a = b = z + 1. 

Proof. We will only prove [Al']. For M ~ {3(1, z) independent on (Vi,V 2 ), the 
characteristic function of the right-hand side of (fTTD with a = b = z + 1 is 


E Jt(AV 1 +BV 2 ) _ Ijg it(l-Jlf/2)Vi , l_ E Jt(V 1 +MV 2 /2) 

2 2 

= 1 f - - / ~ , du + 

2 7„ (l-fl-«/21t0' +1 2 


The last term above is equal to (1 — it) 2 X (1 — it/ 2) 1 /2 due to the well-known 
property of “beta-gamma algebra” that for (A", Y) r\j (3(a,b) ® T a+ fc, one has XY rs./ 
r o . The integral in the second last term can be computed and is equal to (1 — 
it)~ z ( 1 — it/ 2) _1 . Therefore, the characteristic function of AV\ + BV\ is equal to 


(1 — it) 2 (1 — it) 2 1 

2(1 - it/ 2) + 2(1 — it/2) 


( 1 - R )- 2 - 1 


which means that equation ( 1171 ) is satisfied. 


□ 


Remark 5. We note that assertion [Al"] of Theorem 0 extends Theorem 4 in [[50] . 
which proved it in the special case z — 1. The solution to the random equation from 
assertion [Al"] in that case was used to model the distribution of the size of voter 
constituencies in a voter model driven by immigration on the line as discussed in 
that paper. 
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Assertion [A2] under conditions of Theorem [5] gives the limiting distribution of 
an (S 2 , Ri\Ri) version of the CGZ model where the splitting point has density 

( z(l-2u) z -\ 0 < u < 1/2, 

f(u) = l z(2u - l) g -\ 1/2 <u< 1/2, (38) 

[ 0, elsewhere, 

and the larger of the intervals is chosen at each step. In particular, the following 
corollary holds true. 

Corollary 1. Starting with the unit interval, choose a random splitting point Si 
with density (j38j) and select the larger of the intervals [0, Si], [Si, 1]. Continue in¬ 
dependently in this way on the chosen subinterval. Then the limiting location of the 
scheme has distribution /3(z + 1, z + 1). 

Note that Corollary Q] extends the results in |6j and [9], where that claim was 
proved for z — 1. 

2.5 Model 2 

Now consider a more general model where p can have support on an arbitrary subset 
of the whole unit square [0, l] 2 . In this more general model, the nested interval 
scheme, Markov chain { X n } on [0,1] given by (jUJ), and Markov chain {X' n } on M + 
given by (1251) are described respectively in (2A), (2B), and (2C) below, with limiting 
distributions given by assertions [A2]-[A4], respectively. 

(2A) Starting with the unit interval Iq = [0,1], the interval Ii is given by [min(Ai, Si), 
max(Ai,Si)] with ~ p (cf. Section 12711) . For n > 1, the (n + l)st 

interval is chosen from I n using the distribution 

_ f pn if A n > B n , 

^ n+1 \ p' n otherwise, 

shifted and scaled to I 2 (see [M2]), with p Q := p. Recall that if either of the 
conditions [M1]-[M3] on p is satisfied, then the independence property [I] 
holds true. 

(2B) A particle starts from a point x G [0,1] and moves to the point Aix+Bi(l — x). 
The process continues in this way, each new step being independent of the past. 

(2C) A particle starts from a point x' G K + and moves to the point A\x' + B{V\ 
for V\ ~ T a , for some a > 0. The process continues in this way, each new step 
being independent of the past. 

First we will consider a special case of Model 2 discussed in the pioneering work 
of DeGroot and Rao [8j. The result of their Example 6 coincides with assertion [A3] 
in the following theorem. 
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Theorem 6. If (A, B) ~ (3(w,y) <8) /3(y,z) for some w,y,z > 0, then [Al']-[A4] 
hold for (a, b ) = (w + y,y + z). 

Proof of Theorem 0 We will only prove [Al 7 ]. Recall the well-known characteristic 
property of the gamma distribution (see [21]): for (V^ V 2 ) ~ r a ®r 6 , the sum V± + V 2 
is independent of V\/(V\ + Vf). ft follows that 

vrh^+v =u, (39) 

where (Vi, V 2 ) is an independent copy of (Vi, V 2 ). 

By (1391) we have that A Vi ~ IA and BV 2 ~ T y , and so by independence, the 
sum A Vi + BV 2 has the same distribution as Vi ~ T w + y , as required. □ 

In particular, if w = z, then [M3] holds, and so assertion [A2] under the con¬ 
ditions of Theorem [6] gives the limiting location of the (G\, Rj\Di) nested interval 
scheme described in the following corollary. 

Corollary 2. Starting with the unit interval, choose two independent points A\ ~ 
f3(z,y) and B\ ~ (3(y,z ) in [0,1], and select the interval with these end points. 
Continue independently in this way on the chosen subinterval. Then the limiting 
location of the scheme has distribution f3(y + z,y + z). 

Corollary 3. If A ~ /3(w,y) and B = 1, then [A1']-[A4] hold for (a, b) = (w+y, y). 
Proof. To prove [Al 7 ], set B — 1 in the proof of Theorem O □ 

Remark 6. The interval-valued process (2A) generated by fi from Theorem [3] is 
the following (S 2 , R g \Dq) interval splitting model that apparently has not appeared 
previously in the literature. The random rule represented by R G prescribes to choose 
the splitting point S n with the shifted and scaled distribution /3(w, y) (resp. /3(y, w )) 
on I n _i if n is ood (resp. even). The deterministic rule represented by D G prescribes 
to select the right (resp. left) interval if n is odd (resp. even). Then, by assertion 
[A2] under the conditions of Theorem [31 the limiting location of this scheme has 
law 0(w + y,y). 

Remark 7. Now consider the process X := { X. n } n > 0 defined in terms of the forward 
iteration ([2]) as follows. If n is even set X n := X n , otherwise set X n := 1 — X n . 
Then X is described as follows. A particle starts from a fixed Ao = x E [0,1], and 
moves in alternating directions at each step. The particle first moves left to X\ with 
the shifted and scaled distribution /3(y,w ) on [0, X 0 \. In the next step, the particle 
moves right to X 2 with the shifted and scaled distribution /3(ui,y ) on [A 0 ,l], and 
then continues this cycle independently of the past. We have from assertion [A3] 
under the conditions of Theorem [3] that the distribution of X 2n converges weakly to 
/3(w + y,y ) and that of X 2n+ i converges weakly to /3(w,w + y) as n —> oo. Although 
this result extends Theorem 1 in [38] where the special case when w = y = 1 was first 
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proved, it is not a new one. It is a special case of a more general model considered 
on p.2 in [ 21] . 

It turns out we can extend Corollary [3] when w = y = 1 as we will show in our 
next theorem. 


Theorem 7. If (A, B) ~ U[p, 1] (8) Be(l — p), p G (0,1), then [ Al']-[A4] hold for 

M) = ( 2,1). 

Proof. We will only prove [Al']. The characteristic function for the first term of the 
right hand side of (TT7T) is given by 


Ee itAVl 


i r 1 du 
1 ~pJ P (1 - itu) 2 


1 

(1 — it)( 1 — itp ) ’ 


Similarly, for the second term we have that 


Ee <«BV 3= + ijlP = l^P 
1 — it 1 — it 


(40) 


Therefore, by independence, the characteristic function of the right hand side of 
(1T71) is (1 — it)~ 2 . Assertion [Al'] holds and the theorem is proved. □ 

Our next result is a simple consequence of the following lemma. 


Lemma 1. Suppose (A,B, V\,V 2 ) ~ I/[0, y/(w + y)\ ® f3(w, y) ® r 2 <g> P w+y+ i for 
w,y > 0. Then 

AVi + BV 2 ~ iW • (41) 


Proof. The characteristic function for the first term of the left hand side of (14T]) is 
given by 

¥p itAv 1 w + y f v/{w+y) du / ity y 1 

y Jo (1 -itu ) 2 v w + y) ' 

The characteristic function of the term BV 2 is 


Ee itBV2 


B(w,y) 1 



u w ~\i -uy- 1 
(l — itu) w+ y +l 


du = (1 


if )-™- 1 


Jty_\ 

w + y) ' 


and therefore by independence, we obtain that the characteristic function of the left 
hand side of (1411) is (1 — it)~ w ~ 1 as required. □ 

Corollary 4. If (A, B) ~ U[0,y/(y + l)]®f}(l,y), then [A1']-[A4] hold for (a, b) = 
(2, y + 2). 

Proof. Assertion [Al'] follows by setting w — 1 in Lemma HJ □ 

Theorem 8. If (A, B ) = (A, A'B') for ( A ', B') ~ /3(w + y, y)<8>/3(y, w), then [Al']- 
[A4] hold for a = b = w + y. 
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Proof. We will only prove [Al"]. For (Vi, V 2 , V 3 ) ~ T y <g) T w <g) T y , one has that 


(A', B') = 


hi + V 2 hi \ 
hi + hi + hV hi + lij ' 


Therefore it suffices to show that 

d V 2 X + hj 
hi + hi + hi’ 

which by Lemma 1 in [25] is equivalent to showing that 


X(Vi + hi + hi) = XV 2 + hi- 


The last relation follows from Theorem 2 (B) in |TT], as required. 


□ 


Remark 8. Clearly A > B a.s. under the conditions of Theorem [SJ and so setting 

(C, D) := (1 - A + B, B/{ 1 - A + B)), 

it follows that X ~ f3{w + y,w + y) satisfies perpetuity equation fl2Fj) for (C, D) ~ 
/3(2 y,w)®p{y,y). 

Remark 9. Note that assertion [A2] under the conditions of Theorem [8] gives the 
limiting distribution of the following (G\, Rj\Dj) nested interval scheme that has 
apparently not appeared previously in the literature. The random rule represented 
by Rj prescribes to choose the interval with right end point A' x ~ j3(w+y, y ), and left 
end point with distribution /3(y,w) shifted and scaled to [0, A'f\. The deterministic 
rule represented by Di prescribes to select Ii as the interval generated from the first 
stage. We continue independently in the same way on Ii ad infinitum. 

In the special case where w and y are integers, the above nested interval scheme 
corresponds to choosing an interval with end points taken from the uniform order 
statistics C/p), C/( 2 ), • • •, 17(d) as follows. For k = 1, 2,..., |_d/2j, if one chooses Ii = 
[U(k), U[d-k+ 1)] and continues in the same way on Ii independently of the past, then 
the limiting location of this scheme has distribution (3(d — k + l,d — k + 1). This 
is an extension of Kennedy’s scheme in HZ] where the above claim was proved for 
k — 1 . 
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